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AN EXPLICIT SOLUTION FOR THE LAGRANGE MULTIPLIERS 
ON THE SINGULAR SUBARC OF AN OPTIMAL TRAJECTORY 

SUMMARY 


The problem of optimal thrust programming of a single stage rocket- 
powered vehicle to achieve maximum payoff at specified end conditions is dis- 
cussed. The vehicle is assumed to operate in vacuo in an inverse square 
gravitational field. The Lagrange multipliers are analytically determined for 
the so-called singular arc and the mass flow rate to maintain a singular arc 
condition is derived. These Lagrange multipliers are directly related to the 
steering program by algebraic equations. 


INTRODUCTION 


An early problem posed by rocket technology to the mathematical disci- 
plines was the determination of the rocket trajectory which maximized the 
difference between the initial and final values of some function of the trajectory 
end points. This problem, first treated by the calculus of variations, has been 
the subject of numerous research papers. 

A large amount of work was done in attempts to optimize the flight pro- 
file as a function of time under the assumption that the magnitude of the thrust 
was constant. That is, the optimization was concerned only with determining the 
direction of thrust as a function of time, but not the magnitude of thrust as a 
function of time. 

As the theoretical developments became more sophisticated, the case of 
variable thrust magnitude was also treated. In this case, the thrust is assumed 
to be bounded between upper and lower limits; the lower limit may be zero 
(a coasting arc) but the upper limit is usually assumed to be finite so that the 
possibility of impulsive velocity change is excluded. 

Concurrently with these efforts, work progressed in directions which may 
be regarded as perturbational effects of the basic problem. Such effects include 
aerodynamic drag, non-spherical gravitational fields, etc. These effects are 
not part of the present study, and in this report the vehicle will be assumed to 
operate in vacuo over a spherical earth. 


It is well known that if the formulation of the optimal trajectory problem 
allows variable thrusting maneuvers to occur, three types of subarcs, which are 
categorized by thrust level, may occur in the overall flight profile. The first 
possibility is arcs which have constant (non-zero) thrust levels over a non- 
zero portion of the trajectory. The second possible type of arc is the zero 
thrust arc. This possibility can be easily justified on an intuitive basis by 
considering an orbit transfer maneuver with a high thrust vehicle between widely 
separated orbits. 


The third type of arc is the so-called singular subarc of an optimal tra- 
jectory. This arc occurs when the portion of the Hamiltonian which yields the 
direction or magnitude of thrust as a function of time (via Euler's equations) 
vanishes identically. In this case, the thrust level is not immediately available 
from the Hamiltonian and further work must be done to determine the magnitude 
of thrust as a function of time. 


As long as the problem is formulated in the calculus of variations frame- 
work, Lagrange multipliers or equivalent auxiliary variables must be intro- 
duced to determine the direction of thrust as a function of time. Much of the 
difficulty in the numerical calculation of optimal trajectories arises from the 
fact that the initial values of these Lagrange multipliers are not known and the 
resulting split boundary value problem requires repeated integrations of the 
equations which result from the theory. Any technique which solves analytically 
for the multipliers (hence, thrust direction) as a function of the state variables, 
time, and constants of integration would greatly reduce the labor and time 
involved in numerical isolations of optimal trajectories. The above classifica- 
tion of the three types of arcs which may occur along an optimal trajectory is 
particularly well suited to statements about which cases allow such an analytic 
solution of the Lagrange multipliers. The case of zero thrust is simplest, and 
yielded such a solution before the other types. This solution has been known for 
several years [1, 2].* 


*W. E. Miner appears to be the first to have obtained this determination. The 
results of his work are reported in Aero ballistics Internal Note 20-63, Trans- 
formation of the A-Vector and Closed Form Solution of the A-Vector on Coast 
Arcs, George C. Marshall Space Flight Center, Huntsville, Alabama, 1963. 
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The case of constant (non- zero) thrust has not, as yet, been completely 
solved In the sense that is referred to here. A solution for three of six Lagrange 
multipliers in terms of state variables and constants of integration has been 
achieved through the use of a vectorial integral. 

It is the third case that is treated in this paper. It is shown that for 
the case of the singular arc of an optimal trajectory, a complete analytic 
solution for the steering program is possible. The solution is demonstrated in 
the form of three simultaneous algebraic equations. Subsidiary equations then 
yield the additional Lagrange multipliers in terms of the solutions to these 
equations. Questions concerning the conditions under which these arcs actually 
do occur as segments of the trajectory are not discussed. For such con- 
siderations, the reader is referred to References 3, 4 and 5. 


The analytic solution presented for the steering program is applicable to 
numerical trajectory calculations. Using these algebraic equations, it is not 
necessary to integrate the differential equations which normally are used to 
calculate the Lagrange multipliers over that portion of the trajectory which 
fulfills the definition of a singular subarc. 


FORMULATION OF THE PROBLEM 


We wish to determine the extremum value of the difference between the 
initial and final values of some function of the end points, G. For example, the 
fuel expenditure is minimized if we choose G as the difference between the 
initial mass, m Q , and the final mass, m f , i. e. , 


G = m - m . . 
o f 


In order to specify the constraints to which the rocket vehicle is subjected, we 
write the rocket equation as 


V 


T_ 

m 



( 1 ) 
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where 


V = velocity 
T = thrust 
m = mass 

u,p = unit vector in the thrust direction 
g = acceleration of gravity 

The thrust, T, is assumed to satisfy the inequality 

T — T — T 
min max 

The velocity, V, is related to the position vector, r , by 

F=V . (2) 

Finally, the thrust, T, and mass flow rate, me, are related by 

T = -me . (3) 

Using the standard techniques of variational calculus, we now introduce 
vector Lagrange multipliers \ and v and a scalar multiplier cr to form the 
Lagrange fundamental function, F, as 

F = A - (V - T/mu x + g) +F*(F- V) + <r(ih + T/c) . (4) 
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The problem we now wish to examine is to determine the extremal values of 


G* = G + j F dt . 


To initiate this study we first rewrite F in the more tractable scalar notation 
by introducing 


r* = (x, y, z) 

(5) 

ii 

< 

< 

< 

(6) 

X, 

<< 

N 


A 

u^, = (cos y cos 6, siny cos 6, sin o) 

(7) 

g = <S X , gy, S z ) 

(8) 


where y and 6 are control angles (i. e. , direction angles of the vector u^,) 

as shown in the figure below and subscripts indicate components, not partial 
derivatives. 


z 



x 
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in 



We may now rewrite F as 


F = A^V - T/m cos 6 cos y + g ) + A 2 (V “ T/m cos 5 siny + g ) 
x a y y 


+ ~ T/m sin 6 + g ) + iq(x - V ) + v 2 (y - V ) 

^ z x y 


+ v 3 ( z - V ) + <t( T/c + m) 
z 


( 9 ) 


The Euler-Lagrange equations, a necessary condition for an extremal 
path, are 


_d_ / 9F \ _ 9F 
dt \ ByJ 8y. 


( 10 ) 


where the variables y. are V , V , V , x, y, z, y, <5, T, and m. 

i x y z 


we have 


Applying the Euler-Lagrange equation to the variables V , V , and V 

x y z 


Aj + Vi = 0 

(ID 

X 2 + v 2 = 0 

(12) 

\ 3 + v 3 = 0 

(13) 


Defining the ordered set (\ lf A 2 , ^ 3 ) as the vector A? and (iq, v 2 , v 3 ) as the 
vector P'T f we can write equations (11), (12), and (13) as 


A + v = 0 . 


(14) 
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Turning our attention to the variables (x, y, z) for the Euler-Lagrange 
equations, we find 


ag. 


ag. 


ag„ 


, X , “z 

Vi 1 8x 8x 3 8x 


= 0 


ag. 


ag. 


ag„ 


Vo ~ ~ ^2 ~ ^-3 ~T^ L = 0 

2 1 8y ^ 8y 3 8y 


ag x 9g 8g 

h - *i - ^2 , - M it- = o 

3 1 8z 2 8z 3 8z 


These may be conveniently rewritten in matrix notation as 



Using our previous notation for the vectors X and v , and abbreviating the 
matrix by A* we have 

v - A* A = 0 . (15) 

The Euler-Lagrange equation for variable y yields 


T/m(Ai cos 6 siny - \ 2 cos 6 cos y) = 0 . 
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For a powered arc, T/m is not equal to zero. If we assume that 




we obtain 


tany = A 2 /^i 


( 16 ) 


Similarly, for variable <5 


T/m(A 1 sin <5 cos y + X 2 s i n $ siny - X 3 cos 6) = 0 


Again, for T/m # 0, we have 


tan 6 = 




A* cos y + \ 2 sin y 


Eliminating y from this equation by use of equation (16) we obtain 


tan 6 = — — 

+ \ 2 2 


(17) 


It may now be seen that the variables y and 6 have been related to the Lagrange 
multipliers via equations (16) and (17) . 

The Euler- Lagrange equation cannot be applied directly to variable T, 
since T enters linearly. Writing the equation for F as 


F = MV x + g x ) + MV y + g y ) + ^3(V z + g z ) + ^i(x - v x ) 

+ v 2 (y + V ) + v 3 (z - V ) + am + T(a/ c TA/m) =0 

y z 
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(where A has been written for n/a 2 + A 2 + A3 2 ) we see that F can become 
independent of T for certain situations; namely, if the coefficient of T vanishes 


ct/c =f A/m = 0 


If the coefficient of T vanishes, the value of T becomes indeterminate and we 
may write (at the moment) only that T is bounded within the physical limits of 
the problem, i. e. , 

T — T — T 
min max 


The sign ambiguity in the expression cr/c =f A/m =0 can be resolved by 
the Weierstrass test. Originally 


tan y = A 2 /Ai 


( 16 ) 


so that 


cos y = 


±Ai_ 

N/A! 2 + A 2 2 


and 


sin y 


= ± 


^2 

A 2 + A| 


with 


sign (cos y) = sign (sin y) . 


Then 


tan 6 = — 

'v/A 2 + A 2 2 
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giving 


sign (tan 6) = sign (cosy) 

Then 

COS 6 = ± J + + 

and 

sin 6 = ± 1 ~ ^ 

Xi + a 2 2 + a| 

with 

sign (cos 6) = sign (sin 6) 

In order to pursue the question of sign choice further, we may write F as 
F = Aj(V x + g x ) + A 2 (V y + g y ) + A 3 (V z + g z ) + Ui(k - V x ) 

+ 1^2 (y - V ) + 1^3 (z - V z ) + am 

+ T[a/m - i/m(Ai cos 6 cos Y + cos sin Y + A 3 sin 6) ] . 

The Weierstrass condition now requires that for a minimum the function 
E = (t, y, y, y* , y*) = F(t, y* f y*) - F(t, y, y) 
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be positive for all sets y v sufficiently near y and for all sets y* . The last 
member of this equation vanishes under consideration of the strong variations 
of the controls y and 6 since y and 6 do not appear. Thus 


E = F(y* y*) - F(y, y) 



For y^ # y and y^ =£ 6 we have 



so that 


T/m(\ 1 cos 6 cos y + X 2 cos <5 siny + \ 3 sin 6) 


s T/m(Xi cos 6 V cos y* + \ 2 cos 6* sin y *+ \ 3 sin 6*) . 


Choosing 


<5 = 6*, y = y* or 


y = r* 


+ t r 


we have 


cos 6^ cos y + X 2 sin y) — 0 . 


Then 


cos 6 


<s/ Xi + X| ) s 0 


il 



iiiiiiHininiiiiini! i 


or 

sign 'n/Aj + *2 = sign cos <5 • 
Choosing 

y* = y , 6 * = <5 or 6 * = 6 + ir 


gives 


cos 6 cos y + A 2 cos <5 sin T + A 3 sin 6^0 


or 

± A' + ^2 + *3 - 0 


then 



Our function F now becomes 

F = \i(V x + g x ) + A 2 (V y + g^) + A 3 (V z + g z ) 

+ Viix - V ) + y 2 (y - V ) + v 3 (z - V ) + am + T(a/c - A/m) = 0 
x y z 

and the condition for the singular arc is 

a/c - A/m = 0 . (18) 
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The remaining sign ambiguity, ± / A 1 2 +A 2 2 , corresponds to the physical choice of 
initial firing direction. This cannot be expected to be a result of mathematics. 
Further discussions will be restricted to the case of the singular arc. 

Our final variable, m, gives the equation 


b - T/m 2 (Ai cos 6 cos y + As cos 6 sin y + A3 sin 6) = 0 . 
Again, expressing y and 6 in terms of Aj, \ 2 , and A3 we have 


b ¥ T/m 2 A = 0 


which becomes - by our discussion of sign choice 


b - T/m 2 A = 0 


(19) 


Before proceeding to the solution for the multipliers, we note that from 
equation ( 14) 


A = 


so that equation (15) may be written 


A =A*A = 0 


( 20 ) 


The matrix A ;,; is too general for our purposes. We specify a Keplerian 

field 


g = - A* 


r_ 

r 3 
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where p is the gravitational parameter MG, i. e. , the mass of the attracting 
planet multiplied by the universal gravitational constant. From this definition 



3x. 8x. 

1 1 


V- 


(x 2 + y 2 



0.. OX.X. 

- ^ (x 2 + y 2 + z 2 ) 3 / 2 " ^ (x 2 + y 2 + z 2 ) 5 / 2 


6.. 3x.x. 

= M ^ M (i.j = 1,2,3) 


where 


x. = x or y or z 

r has been written for ^ x 2 + y 2 + z 2 


6.. is the Kronecker delta, 
ij 


rv ^ )r 

We can now write the expression A V X under the Keplerian field assump- 
tion, A \ , as 
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3xyA? 

r 5 



= M 


3*y^i | ^2_ 
r 5 + r 3 



3yzAa 



3yzX, 
r 5 + 


Aa_ 

r*3 



ft 


/^T (xAi + yA 2 + zA 3 )\ 

r 3 

, 

-3// 

(^l + y^-2 + z ^3> 

w 


(xAi + yA 2 + zA 3 ) j 


X „ A • r" -A 

= V " 


Equation (20) may now be written as 


A + M 



= 0 


We may also rewrite the basic equation of motion, equation (1) , into a 
form not involving the angles y and 6 by use of equations (7) , (16) , and (17) . 
Thus 


T A , r 
r - m A + 


= 0 
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The principal results of this section may now be summarized as 


rr T X 
r - — — + 
m X 



( 21 ) 


X + 





= 0 


( 22 ) 


a - 



(23) 


£ 

c m 


0 


Equation (24) specifies that we remain on a singular arc. 


(24) 


Before proceeding to the integration of this system of equations, it is 
worthwhile to consider the following. In equation (21) we have three position 
and three velocity components which must be initially specified. Inequations 
(22) and (23) we must specify the initial values of the vectors X, X , and of the 
scalars a and T (or cr). This gives 14 initial conditions. Since the equations 
in the Lagrange multipliers are homogeneous, one initial value is arbitrary. 
Furthermore, equation (24) gives a relationship between X and a which must 
be maintained, thus reducing the number of initial conditions to 12. For a 
solution which yields the Lagrange multipliers as an algebraic relationship, 
we may expect to obtain six independent constants of integration. 


DERIVATION OF INTEGRALS 


From equation (23) we have 
TX 

m<7 - = 0 , (25) 

m 
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and from equation (3) 


mu + = 0 . (26) 

c 

Using equation (24), equation (26) becomes 
TA 

ma + — = 0 . (27) 

m 

Adding equations (25) and (27) we obtain 
* (m ' 7) = 0 


or 


ma = Cj (28) 

where is an integration constant. 

To obtain our second integral, we note that the Euler-Lagrange equations 
are explicitly independent of time. From the general theory of variational 
calculus we know that 

10 

F - Z h 1 * = °» 

i=l y i 

where 

F is given by equation (9) 

y. represents V x » V , V z , x, y, z, y, 6, T, or m. 
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Furthermore, since F is identically equal to zero, the above equation becomes 


10 

-E *i F y = C 2 

i=i y i 


By differentiating and summation we obtain 

\lV + X 2 V + X 3 V + xqx + V 2 y + + o™ = - c 2 

x y z 

• • • • • • 
Substituting V = x, V = y , V = z , v 2 = -X t , v 2 = -A 2 , and i> 3 = -X 3 , this 
becomes x y z 

/ • • • 

Xix + X 2 y + X 3 z - ^x - Xgy - X 3 z + am = - C 2 , 


which may be written as 


X • r - X • r + am = - C 2 


(29) 


Taking the dot product of X with r as given by equation (21) and inserting 
the result into equation (29) we find 


X • r + n X ‘ r T £ " “ ( m X + ar ”) = C 2 ' 


The term in parentheses of this equation may be written 

T / TX \ A . 

— X + am = f — 0 1 nx + am 
m y m ) 


= am + am 


= dt (m<T) = 0 • 
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where equation (23) has been used. Thus, our second integral is 


X 


r + n 



— C2 . 


( 30 ) 


Following PLnes [ 6] , we begin derivation of the third integral by noting 

that 




r) = X • T+’X 


= X • r + 



= X • r + 2 p 


Substituting X • r from equation (30) we have 


_d 

dt 


(A. * r ) = C 2 + M 



(31) 


Taking the dot product of equation (21) with X, the dot product of equation (22) 
with rand subtracting, gives 


— (if 

dt 


r - X • r)= - — X + 3 u 
m 
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so that 




3 dt (A 


-r 7* t*. T 


r - X • r ) + 


3m 
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A, * r 

Inserting p — pj — frona this expression into equation (31) gives 


d 


— , „ i d 

r) =C 2 + Iii 


r - A 


F)+ a 

3m 


or 

~ (2A • r + A • r ) = 3C, + — A . 
dt m 


Now, using the defining equation for T in terms of m, equation (3) , this may 
be written as 


_d 

dt 


(2A • r + A • r ) = 3C 2 - cA — 


(32) 


To perform the integration we must first note that equations (24) and (28) 

yield 

A = -^ (33) 

c 

making A a constant. Equation (32) now integrates to yield 
t*' — *• — ► — *• 

2A • r + A • r = 3C 2 t + Cj H n -pp + C 3 . ( 34) 

To obtain our final integral, we return to equations ( 2 i) and ( 22 ) . From 

these 

V . A x r 
A x r + p g — = 0 

r- vr A x r 

A x r + u 3 — = 0 

r 
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Subtracting gives an immediately integrable equation yielding 


Axr-Xxr = C 4 


where C 4 = (C 4 , C 5 , C 6 ) is a vector integration constant. 

To determine if the constants C 4 , C 5 , and C 6 are independent, we apply 
the Jacobi test. Thus 


C 4 = - \ 3 y - \ 2 z + *3y 


Cg — A 3 x - X 4 z — A 3 x + A 4 z 


c 6 = ^iy - ^2 x - My + a 2 x 


0 . ^ = -^3. ^ = ^2 . ^r = o , ^ = A 3 , ^=-A 2 


ac K • ac^ ac 

ax 3 ’ ay ’ 8 z 


ac s : ac5_ acs _ A ac fi , 

az " _A * * ak " ~ As ’ ay " 0 ’ az " Al 




we may write 


0 9 ^3 9 A 2 > ^ > A 3 , -A 2 

J = Ag , 0 , -Ai , -A 3 9 0 , A^ 


-Ag 9 A^ , 0 , A 2 , -Ai , 0 



A non-vanishing sub-determinant is 
• • 

-A 3 > ^2 » 0 

0 , -A-i , A 3 = — Aj ( Aj> A 3 - A 2 A. 3 ) ^ 0 
Aj , 0 , A 2 

in general, so our integration constants C 4 , C 5 , and C 6 are independent. 

The preceding statement cannot be taken to imply that we may solve for 
^ 2 » a nd A 3 in terms of A , r, r, and C 4 . To see this we write equation (35) as 



Now 

0 z -y 

A = -z 0 x = zxy - zxy = 0 
y -x 0 

• • T 

so that we cannot invert the coefficient matrix of (A*, A^ A 3 ) . This situation 
occurs in similar fashion with the angular momentum integral of the three -body 
problem. 

DETERMINATION OF THE LAGRANGE MULTIPLIERS 


The purpose of this section is to demonstrate three independent equations 
which involve Lagrange multipliers Aj, A 2 , and A 3 algebraically. Ihe first such 
equation, equation (33), has already been derived. The second will be derived 
from equation ( 35) . From this relationship we obtain 

r* • ( t* x A) = r* • C 4 (36) 
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since r • (X x r ) = 0 . 

To obtain our third relationship we first note that from equation (35) we 

obtain 


r x ( X x r) = r x ( C 4 + X x r ) , 
or 

( r • r ) X =(r • \)r + rx (C 4 +Axr) 


Inserting r • X from equation (30) into this expression gives 


( r • r ) X = \ C 2 - m 


X • r 


r + r x ( C 4 + X x r ) 


(37) 


Now taking the cross product of r with equation (35) gives 


(r*r)X=(r-X)r + rx(C 4 +Xxr) . 


Introducing r • X from equation (34) gives 


(r • r)X = - 


T* 1 / “O 

x = - C t in 

2 \ m 


+ 3 C 2 1 + C 3 — X 




+ rx(C 4 +Xxr) 


(38) 


Eliminating X between equation (37) and (38) under the assumption that 
r • i* =£0 gives 


[}(«<■=! 


r) 


+ 3C 2 t+C 3 - A* r) r+rx(C 4 + A x r ) ( r • r) 


■] 


= £ (c 2 - M - r 3 r - )r r +r'x(C " 4 + Xx r) J ( r • r ) . 
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This equation can now be regrouped to give the form 


A • r \ 1 t* v 

Qj -M —^3— ) ~2 {r ' r } 

+ r x (C 4 +A xr)(r * r) - r x ( C 4 + A x r ) ( r • r) = 0 . 
The last two members of this equation can be written as 

r x C 4 ( r* r) -rxC 4 (r- r) + [ r ( r • r) - r ( r • r)]x ( A x r 
= -C 4 x[rx(rxr)] + [rx(r x r ) ] x ( A x r ) 

= [r x (r x r ) ] x (C 4 + A x r ) = [r x (r x r ) ] x (A x r) 

= - (A x r) • (rxr)r + [(A x r) • r](r xr) 

= - (A x r) ' (r x r)r 
where we have used 

(A x r ) • r = 0 . 





111 _ 
C 4 in -^+3 C^t + Q - A * F 


Using the final result, equation (39) may be written as 


r 2 (c 2 - 


A • r \ 1 t*\ 

M — — I" o < r ' r ) 


in ^ +3 C 2 t + 03 ^ 


Since 0, the coefficient must vanish. Introducing A x r back into the 
coefficient yields 


r 


2 





r ) 



m 

n + 3 Cot + 
m l 


C 3 


- A • 



- ( Qj + \ x r”) • (iT x r) 


0 


(40) 


Equations (33), (36), and (40) now define the Lagrange multipliers 

A 2 , and \ 3 in terms of r, r, m, and the integration constants. The 
multiplier cr comes from equation (28) and iq, iq, iq can be obtained from 

equations (11) , (12) , and (13) once X is computed from equation (35) . 


DERIVATION OF THE MASS FLOW RATE EQUATION 


The preceding information is not sufficient for a complete solution of 
the problem for we do not yet know either m or T as a function of the problem 
variables. The necessary information is available, though obscure. 

From equation (24),. which is the requirement that the singular arc be 
followed, we have 


a 

c 


X. Am 

+ 2 

m m 


= 0 


Inserting the values for a and m we find that 


A = 0 

or A is a constant which we have ascertained from equation (33) . Furthermore 
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so that 



From equation (22) we obtain 



The derivative of equation (42) is 



but from equations (41) and (22) 

t * 7f X • r T * — *• 

X • X = 3m-^T~ (A • r) 


( 41 ) 


(42) 


(43) 
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so that equation (43) may be written 


[(f)' 


r‘ II —| + 2 ( r- F)(X- F + 2 X 


A • F)J 


5(A • rV(r ■ F) = 0 


By use of equation (34) , the previous equation may be written 


(f) 


m 


Cl ^ (F* r )+ 2(A • F)^3C 2 t + Cj in -£■ + C 3 J 


- 5(A • r ) 2 (r • r ) = 0 . 


(44) 


The derivative of equation (44) is 


2(r • r ) 2 A 2 + 


i n + 3 C 2 t + C$j | r 2 ^ A • i 


m 


r + C t in — “ + 3C 2 t + C, 
1 m it 


) 


(A ' r ) (r • r) - 5(A * r ) (A- r ) (r • r) + 6 C 2 r 2 (A • r) 


+ [ r 2 A 2 - 5( A • r ) 2 ] I r • r - 


[*• 


k i 
rj 


= (m)vr) [3xV - 5(r ' ^ 2] 


(45) 
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Now c * 0 for a powered arc, and assuming 


X ' r * 0 

we can solve for m, with m given by equation (441, if 

3 A 2 r 2 - 5(T- 7 ) 2 * 0 . (46) 

From equation (42) , as noted by Corben [7] , we may write 

r 5 T- r = /u[X 2 r 2 - 3(X- 7) 2 ] . 

Since the left side of this equation is non-negative, we can write the angle be- 
tween X and r as a and conclude 

i - 3 cos 2 a 2: 0 
or 

^ s cos 2 a . 

Our criterion of equation (46) becomes 
cos 2 a * 3/5 . 

Since this condition is excluded by the previous inequality, we need only insure 
that 

T- * 0 

to solve for m; m is given by equation (44) . 
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CONCLUSIONS 


The general problem of an optimal rocket trajectory in vacuo in an 
inverse square gravitational field is composed of three possible types of seg- 
ments. The first of these is the coasting arc, and, for this case, a complete 
determination of the Lagrange multipliers as algebraic functions of the state 
variables, time, and integration constants has been obtained [1, 2]. 

The second case, sub-arcs of constant thrust, is most important from 
the practical standpoint and has not, as yet, yielded an analytical solution. 

The third case is sub-arcs of optimal intermediate thrust (singular 
sub-arcs) and this area has been treated in this report. It has been shown that 
integrals which exist in the literature, equations (28), (30), (34), and (35), are 
sufficient to explicitly yield the set of Lagrange multipliers which govern the 
guidance law along an optimal trajectory. The simultaneous algebraic equations 
governing the Lagrange multipliers associated with dynamical constraints are 
equations (33) , (36) , and (40) . Once this set of equations has been solved, 
equations ( 11) , ( 12) , ( 13) , and ( 28) then yield the additional multipliers which 
were associated with kinematical constraints in the original formulation. Equa- 
tions (16) and (17) relate the Lagrange multipliers to the thrust direction angles. 
Finally, equation (45) supplies necessary information about the mass flow rate 
which is required to maintain the singular arc condition. Using these equations 
it is possible to compute an optimal singular subarc without encountering the 
difficulties associated with the usual split boundary value problem. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, February 20, 1967 
933-50-01-00-62 
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